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ABSTRACT 

The goal of searching back in cosmic time to find star formation during the epoch of reionization will soon 
be within reach. We assess the detectability of high-redshift galaxies by combining cosmological hydrody- 
namic simulations of galaxy formation, stellar evolution models appropriate for the first generations of stars, 
and estimates of the efficiency for Lya to escape from forming galaxies into the intergalactic medium. Our 
simulated observations show that Lya emission at z ~ 8 may be observable in the near-infrared with 8-meter 
class telescopes and present-day technology. Not only is the detection of early star-forming objects vital to 
understanding the underlying cause of the reionization of the universe, but the timely discovery of a z > 7 
star-forming population — or even an interesting upper limit on the emergent flux from these objects — will 
have implications for the design of the next generation of ground- and space-based facilities. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Recent observations have significantly advanced our un- 
derstanding of the reionization of the universe. Estimates 
of the Thomson optical depth from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) suggest t hat the reion- 
ization epoch began at red shifts 14 < z < 20 jBennett et alJ 
120031 ISpergel et alJ |2003) and Gunn-Peterson troughs in 
distant quasars indicate that reionization ended at z ~ 6 
jBecker et alJl2001t iFan et alJ 120021) . However, the known 
populati on of quasars ca nnot produce sufficient ionizing ra- 
diation iFan et alJ 120 01). so star-forming galaxies are most 
likely the dominant source of ionizing photons at redshifts 
z > 7 . a view supported by theoretical st udies of reionization 
(e.g.. lTinslevlll973r lSokasian et al. 2003a). The direct detec- 
tion of these objects lies at the next frontier in the study of the 
evolution of the early universe. 

Discoveries of star-forming objects are progressing to 
ever higher redshifts, with successfu l results from both 
broad- band color selection (e.g.. lYan. Windhorst. & Cohen] 
1 2003 1 : IStanwav Bunker. & McMa honl l2003t iBunker et alJ 
120031 [Dickinso n et alJl2003l) and Lya emission li ne searches 
(e.g. iRhoadsetalJ I2003t iMaier et alJ 120031 iCubv et alJ 
2003). Although the number counts of /-band "dropouts" 
(at 5.5 < z < 6.5) are still heavily debated, it is clear that a 
substantial population of galaxies with detectable continuum 
exists at these redshifts. The next step — discovering a popu- 
lation of objects before the end of reionization — lies behind 
a technological barrier: most of the light from z > 7 objects 
falls in the near-infrared, where backgrounds provide signif- 
icant challenges from the ground. High spectral resolution 
observations provide some relief; dispersing the light elimi- 
nates much of the backgr ound noise and allows the detection 
of lin e em ission (e.g.. | T hompson. Diorgo vski. & Traugerl 
19951 iThomnson. Mannucci. & Beckwithl 119961: 
Panaris. Stiavelli. & Falll2003l~ 

At first glance, weak star formation and a neutral intergalac- 
tic medium (IGM) would appear to critically hamper the de- 
tection of Lya at high redshift. However, a top-heavy initial 
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stellar mass function (IMF) and l ow metallicity may boost the 
Lya flux of the earliest stars (e.g., Bromm. Kudritzki. & Loeb 
1200 It ISchaererll2003l) . In addition, this emission ma y pene- 
trate even a neutral IGM with reasonable efficiency ( Haiman 
2002: :Santosl2003l) . 

Taken together, the arguments above strongly motivate the 
search for Lya emission at redshifts z > 7 using narrow-band 
imaging in the near-infrared. Here, we describe a search for 
Lya in a narrow /-band "window" in the sky background, 
corresponding to z ~ 8. 

2. STAR-FORMING OBJECTS AT Z sa 8 

2.1. Hydrodynamical Simulations of Early Star Formation 

We examine th e star formation rates in the cosmologi- 
cal simulations of Springel & Hernquist (2003a). Employ- 
ing a novel treatment of smoothed particle hydrodynam- 
ics and a multiphase description of star-formin g gas to in- 
corporate fee d back and stellar wind s ( Springel & Hernquist 
120021 12003bl) . Springel & Hernquis] J2003al) obtained a nu- 
merically converged prediction for the cosmic star formation 
rate (SFR) as a function of redshift that agrees with the ob- 
served star formation history at low redshift. 5 Here, we fo- 
cus on the 10 h~ l comoving Mpc Q5 simulation of a WMAP- 
concordant cosmology with 2 x 324 3 particles. The more lu- 
minous sources in the simulations have star formation rates 
in the range of 10~ 2 < M < 1 M /yr. With typical dynami- 
cal masses of 10 8 - 10 9 M at z = 8, their intrinsic circular 
velocities are generally 30 -50 km s _1 . The simulations pre- 
dict an early rise in star formation activity at z > 5; at z = 8, 
the star formation rate density is within a factor of ^2 of its 
peak, relatively independent of the details of the physics of 
star formation and feedback (Hernquist & Springel 2003). 

2.2. Lya Emission and Absorption: the IMF, the IGM, and 

the ISM 

The amount of Lya that we can observe from young star- 
forming objects depends on the IMF and metallicity in early 
galaxies, the amount of attenuation within the galaxy itself, 
and the amount of scattering in the intergalactic medium: 

FLya.obs = e" Tdust ( 1 - f esc ) flGM F Lyctjem ( 1 ) 

5 We note an error in figure 1 2 of SDrinael & Hernauist 1200.3 all in which 
the observational estimates of the SFR were plotted too high by a factor of 
h = 1.4. When corrected, the observed points are in better agreement with 
their theoretical estimates. 
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where F^ ya b s is the observable Lya flux per unit of star for- 
mation, Tdust is the optical depth due to dust in the galaxy, 
f esc is the fraction of ionizing photons that escape the star- 
forming galaxy and thus create no Lya photons at the source, 
and fiGM is the total fraction of Lya photons that penetrate 
the intergalactic medium and reach us. FLy Q , e m is the maxi- 
mum, ionization-bounded flux of Lya photons produced per 
unit of star formation; it depends on the IMF and metallicity. 
Recent models and the interpretation of high-redshift galaxies 
provide some constraints on these quantities. 

The IMF: Theoretical models predict a top-heavy 

IMF for the first generation of star s form ed in 

metal-free gas JBromm. Coppi. & Larsonl IT9991 |2002; 
Abel, Bryan, & Norman 2000). These stars can radiate more 
than an order of magnitude more ionizing photons per unit 
solar mass than a Salpet er IMF, resulting in much stronger 
Lya and HeI I(A1640) JBromm. Kudritzki. & Loebl 1200 It 
ISchaererl 120031) . However, only the initial generations of 
stars will be formed from completely pristine gas; thus, one 
major source of uncertainty is whether luminous bursts of 
star formation later in the process (z = 8) will be too enriched 
to exhibit extremely top-heavy IMFs. 

The ISM: The factor e^** (1 — fesc) is the fraction of emit- 
ted ionizing photons available for conversion to Lya; it in- 
cludes the absorption of ionizing and Lya photons from dust 
near the newly forming stars. There are, at present, few di- 
rect observatio nal cons t raints on this factor at any redshift. 
At z ~ 3, Petti ni et alJ ( 119981) estimate e" Tdust ~ 0.16-0.4. 
But the stellar populations at even higher redshift are prob- 
ably younger, less chemi cally evolved, and hence less dusty . 
For the escape fraction, Steidel, Pettini, & Adelberger (2001) 
combine z ~ 3 Lyman-break galaxy spectra to estimate an av- 
erage f esc > 0.07-0.1. 

The IGM: The discovery of Gunn-Peterson troughs in the 
spectra of quasars at z > 6 indicate that the average neutral 
fractio n of the intergalactic medium exc eeds xhi ~ 10~ 2 at 
z > 6 llBecker et alJ200lUFan et all2002D . At face value, this 
result suggests that Lya emission from beyond this r edshift 
will not penetrate the IGM. However, Haiman (2002) shows 
that a substantial fraction of the emitted Lya photons escape 
scattering by the IGM when a galaxy ionizes a local bubble 
in its immediate surroundings. Extending this idea using dy- 
namical models of the IGM and galactic winds, Santos (2003) 
derives values of fiGM ranging from < 0.002 to ~ 1. This 
broad range illustrates the substantial uncertainties involved 
in estimating the expected Lya flux. In general, penetration 
through the IGM improves for galaxies with lower redshifts, 
higher star-formation rates, and older (longer-duration) bursts 
of star formation, higher f esc , outflows in which the centroid of 
Lya is shifted redward, and with higher ionizing backgrounds 
in the universe. Because unabsorbed ionizing photons either 
create Lya in the source or escape to ionize the IGM, but 
not both, the total emergent Lya decreases for both ver y high 
and ve ry low escape fractions. The fiducial models of Santos 
(2003) show a broad peak in detectable Lya at f esc ~ 0. 1-0.8. 

The strongest constraints on the creation and transmission 
of Lya come from the cosmic microwave background. The 
WMAP observation of th e Thompson op tical depth to reion- 
ization, r e = 0. 17 ± 0.4 ( Kogu t et al.l2003l) . constrains the pa- 
rameters in Equation n By modeling the reionization pro- 
cess, Cen (2003a,b) and Sokasian et al. (2003b) find that a 
large Thomson optical depth requires most or all of the fol- 
lowing at high redshift: a top-heavy IMF of extremely low- 
metallicity stars, a high star formation efficiency in low-mass 
halos, f esc > 0.3, a positive tilt in the matter power spectrum, 
or an additional source of ionizing photons. Except for es- 
cape fractions extremely close to unity, these conditions are 
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FIG. 1. — "Windows" in the /-band night sky spectrum. The black line 
indicates the transmission of the night sky (scale on right). We plot the night 
sky spectrum at two resolutions (R = 1000, red line; R — 300, blue line) to 
indicate regions below 1/4 of the mean background (blue shading for moder- 
ate resolution, red shading for high resolution). At moderate resolution, only 
a small number of windows are available. 



all conducive to observing z > 7 galaxies in Lya. 

Although much uncertainty surrounds the search for Lya in 
z > 7 galaxies, we argue that the distinct possibility of suc- 
cess motivates an attempt with present-day technology. As an 
illustration, we consider two specific scenarios that are both 
detectable in Lya by design. They differ by a factor of ^3.4 
in the flux of observable Lya photons. They are: 

1. Plausible: An IMF rich in massive stars with a 
Salpeter slope, containing only stars in the range 50 - 
500 M , with a metallicity of Z = 0, an escape 
fraction of f esc = 0.1, T(j ust = 0, and fiGM = 0.25, 
which is appropriate for a range of IGM wind mod- 
els, yielding F , v „ nh . = 6.4 x 10 42 erg s" 1 / [M ^ /yr] of 
detectable Lya (JSchaerer 2003; Santos 200|). 

2. Optimistic: An extremely top-heavy IMF com- 
posed of only 300-1000 M R ze ro-metallicity stars 
( Bromm, Kudritzki, & Loeb 2001), Td ust = 0, f esc = 
0.35 to maximi ze the observa ble Lya emission based 
on the fiducial ISantosI |2003) models, and fiGM = 1, 
appropriate for an IGM model with some ionizing 
background in which strong galactic winds clear the 
Stromgren sph eres around galaxies and shift Lya to- 
ward the red (Santos 2003). This combination yields 
FLya.obs = 2.1 x 10 43 erg s _1 /[M /yr] of detectable 
Lya. 

Both scenarios make assumptions that are relatively favor- 
able for detecting Lya emission, such as the presence of a 
top-heavy IMF in high-redshift galaxies. If regions of intense 
star formation are sufficiently chemically enriched by z ~ 8 to 
have a Salpeter IMF from 1-100 M Q , FL ya ,em is reduced by 
a factor of 4- 12, depending on metallicity. With 6.3 x 10 8 
years of cosmic time available for star formation before z = 8, 
delaying chemical enrichment sufficiently long may require 
suppressing star formation at early times, as would be the 
case if the power spec trum were reduced on small scales (e.g., 
Yoshida et al. 2003a b). In addition, inhomogeneous structure 
formation may lead to pocket s of low-metalhcity star f orma- 
tion e ven at late epochs ( Scanna pieco. Schneider. & Ferraral 
2003). Observationally, rest-frame equivalent widths of Lya 
as late as z ~ 5.7 support the hypothe sis that the IM F may 
be top-heavy at late epochs (e.g., Rhoads et al. 2003). More 
indirect arguments in favor of massive and/or low-metallicity 
early stars come from comparisons of the star formation his- 
tory of the universe to the requirements for reionization (e.g., 
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FIG. 2. — A realization of a 32-hour 2' x 2' observation of the "optimistic" 
scenario with an 8-meter telescope. We simulate detection of the objects 
in the Sorinael & Hernquist 1 2003a) Q5 10 IT 1 comoving Mpc model. The 
R = 125 narrow-band filter we assume is actually 44% thicker in redshift 
depth than this simulated observation. There are 56 detectable sources in the 
field. We show the dotted inset enlarged for clarity. 
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FIG. 3. — Lya fluxes in the z = 8.2 27 /-band window. We plot num- 
ber counts per square arcminute in the Snringel & Hernquist 1 2003a) Q5 
101r' comoving Mpc simulation, converting star formation rates to observ- 
able Lya fluxes via the "plausible" and "optimistic" scenarios outlined above. 
The arrows indicate detection limits for 32-hour exposures with 8-meter and 
30-meter telescopes at moderate (R = 125) and high (R = 3000) resolutions. 



Ferguson. Dick inson ., <fe Pap ovich|| 2002) and from fossil ele- 
mental abundances (Beasley et al. 2003). 

3 . IDENTIFYING AND MAPPING HIGH-REDSHIFT SOURCES 

A simple argument shows that for even a moderate Lya 
emitter, detecting Lya from the ground at high spectral res- 
olution is easier than detecting the continuum emission. For 
example, in the J-band window at 1 .122fim (Lya at z=8.2; see 
Sec. 3.1), assuming a Mauna Kea sky background and obser- 
vations in which sky noise dominates, the signal-to-noise ra- 
tio of a suitably located R = 200 observation of Lya becomes 
equal to or better than that of a broad-band continuum obser- 
vation (R = 4) when the rest-frame EW(Lya) em > 27 A. At 
R = 2000, the threshold is EW(Lya) em > 8 A. Observed rest- 
frame Lya equivalen t widths of z = 5.7 galaxies exceed 150 A 
jRhoads et alJl2003l) and starburst models predict equivalent 
wid ths far in excess of 100 A for most young starburst scenar- 
ios (iBromm. Ku dritzki. & Loeb 2001; Schae rerfcOOl . Thus, 
detection in Lya is likely the most powerful way to identify 
and map very high-redshift star formation from the ground. 

Observing spectral features between the strong atmospheric 
lines helps to minimize the background in the near-infrared. 
Although there are existing spectrographs with high enough 
resolution to observe between OH lines, the relatively low 
number density of detectable objects predicted by simulations 
argues for blind searches for Lya emitters in a large, con- 
tiguous field (> few □')■ Because narrow-band filters with 
R > 200 are not widely available, the most effective current 
strategy is to identify the broadest OH-free windows in the 
night sky; Fig. ^ illustrates some windows in the /-band, in- 
cluding a wide window at z = 8.2-8.3 with ^80% (variable) 
atmospheric transmission. 6 

3.1. A Simulated Observation with an 8-meter Telescope 

We assume each star-forming particle in our simulation is 
a source of Lya photo ns, c onverting star formation rate to 
FLy Q ,obs as described in 32.21 For an R = 125 narrow-band fil- 
ter in the z = 8.227 /-band window, 0"35 seeing, a total sys- 

6 We use the Gemini sky spectrum from http://www.gemini.edu 



tem throughput of 20%, and the detector parameters appro- 
priate for a modern detector (the NIRI detector on Gemini), 
we simulate 32-hour observations of Lya for the "optimistic" 
and "plausible" scenarios. Fig. [2] shows a 2' x 2' corner of 
the "optimistic" simulation; the thickness of the simulation 
is 44% less than the thickness of an R = 125 observation, so 
the simulation actually underrepresents the number of sources 
expected in this scenario. 

Allowing for the depth and angular size of the simulation, 
the volume probed by the entire Q5 model is 5. 14 times larger 
than the volume of a single R = 125 narrow-band filter ob- 
servation with a 2' x 2' detector. We analy ze the simulated 
obser vations with the SExtractor software (Bertin & Arnouts 
1996). For the "optimistic" case, we find 176 sources in the 
entire simulation volume to a S/N of > 3 .4. Although the seg- 
ment in Fig. 13 is in a well-populated part of the simulation, 
with 56 detectable sources, the detectable objects do cover 
the volume well in the sense that any 2' x 2' field would in- 
clude observable galaxies. In the "plausible" simulation case, 
we detect 23 sources. The objects are clustered, but a 2' x 2' 
observation would reveal no sources <20% of the time. 

3.2. Identifying Low-redshift Contaminants 

Although detecting galaxies via their Lya emission is 
the primary challenge, separating the sources from possi- 
ble contaminants is equally important. Both continuum and 
emission-line sources will appear in the narrow-band images. 
Contaminants from lower redshifts must be identified with ad- 
ditional data. Deep broad-band images in both the optical and 
near-infrared will serve this purpose. For example, for a 32- 
hour narrow-band observation at R = 125, /-band data that is 
deeper than ^25.9 (AB) should reveal all continuum sources 
in the narrow-band image. Many upcoming imaging surveys 
will exceed this sensitivity. 

Extremely high redshift galaxies may be difficult to detect 
in the c ontinuum, even wit h HST The brightest galaxy known 
at z ~ 6 ( Dickin son et al!20 " Q3"t) would be detectable at z = 8 .2, 
without evolution, in the upcoming HST/NICMOS ultra-deep 
field observations. But at z ~ 8 galaxies should have lower 
star formation rates and less stellar mass. Even if the contin- 
uum flux is detectable, the target objects will be identifiable as 
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z-band dropouts and distinguishable from lower redshift con- 
taminants which will appear in the optical bands. 

The remaining possible contaminants are low-z objects 
with strong emission lines and no detectable continuum. In 
all cases, however, the emission lines must be prohibitively 
strong for the objects to remain undetected in extremely deep 
broadband images. For the scenario described above, if a 
detected emission line in the 1.122 /im window is actually 
[O II](A3727), the z = 2.01 galaxy would have to have an un- 
physical equivalent width of ~800 A or more (rest-frame) to 
remain undetected in a deep F606W observation (29.3 magni- 
tude limit at 3er, corresponding to the HDF-North data). Other 
emission lines place similar requirements on the high-redshift 
objects, including both QSO lines and strong galaxy lines. 
Thus, the sources that appear in the narrow-band image but 
are non-detections in deep optical continuum images are al- 
most certainly z > 7 galaxies. 

4. GOALS FOR PRESENT AND FUTURE FACILITIES 

Many far-reaching science goals motivate devoting an enor- 
mous effort to the study of extremely high-redshift "first- 
light" objects. Some of these goals are: (1) measuring the 
luminosity function of z > 7 galaxies, to compare with the 
number of ionizing photons needed to reionize the universe, 
and to measure the global star formation rate at these red- 
shifts, (2) mapping high-redshift sources and measuring their 
kinematics, to observe their clustering properties and compare 
with models of early galaxy formation, (3) quantifying high-z 
Stromgren spheres in the IGM from Lya line profiles, and (4) 
directly probing the IMF and metallicity of early stars with, 
e.g., He II(A1640) and, eventually, other spectral lines. 

Accomplishing these goals will require a large, comple- 
mentary set of observational facilities. The James Webb Space 
Telescope, with broad-band sensitivities of ^2 nJy, will detect 
early stars in their continuum light. Narrow-band imaging 
will likely remain the strongest tool of ground-based facilities 
in the near-infrared. Fig. |3| shows the expected performance 
of 8-meter and future 30-meter ground-based telescopes. As- 
suming modern detector characteristics, we plot the fluxes 

7 This technology is one goal of the DAZLE spectrograph, under develop- 
ment for the VLT i http://www.aao.gov.au/dazle/ 1. 

8 At, e.g., R = 3000, a long NIRI observation would actually be detector- 
noise dominated; thus, the development of low read noise and low dark cur- 



from star-forming regions at z = 8.227. Clearly, 8-meter class 
telescopes have great near-term potential for advancing the 
study of early star formation. Narrower filters (R > 1000) 
would improve sensitivity substantially, although they would 
require some development in both filter and detector technol- 
ogy. 1 ■■ 8 Looking further into the future, the detailed physical 
study of these objects — an exploration down the luminosity 
function, measurements of Lya line profiles, and the detection 
of fainter spectral emission lines — is probably the role of the 
next generation of large ground-based telescopes. The timely 
discovery of high-z Lya emitters is crucial to guide the design 
of future facilities, since probing the epoch of reionization is 
one of their central scientific goals. 

In this Letter, we motivate the search for high-redshift 
galaxies by presenting plausible scenarios in which the Lya 
emission from z ~ 8 galaxies is observable with 8-meter-class 
telescopes. The discovery of z > 7 objects is a crucial com- 
ponent of the quest to understand the evolution of the early 
universe. They existed at an epoch when the universe was 
still partly neutral, and they are the best candidates for sup- 
plying the radiation that reionized the universe. Although we 
are able to combine many areas of recent research to arrive at 
an educated guess about the properties of z > 7 objects, there 
are still many critical unanswered questions about the nature 
of reionization. With the discovery and study of z > 7 ob- 
jects, many of these questions will be answered. We stand on 
the verge of probing reionization directly. 
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rent detectors is another important component in the search for extremely 
high-redshift objects. 
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